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INITIAL DATA 



Material properties: porosity, water adsorption isotherm, 
pore solution, solid phase content, 
water content, diffusion coefficients, 
water diffusivity, ion adsorption 
isotherm, temperature 

Numerical parameters; number of rime steps 

lengch of time steps, FEM mesh 

Boudary conditions: concentrations, water content, 

electrical potential 
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TRANSPORT STEP 

Resolution of equations (2), (7), (14) with the finite 
element method ^ 
OUTPUT: concentration profile for each ion 
water content profile 
electrical potential profile 
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CHEMICAL REACTION STEP 

Resolution of the chemical equilibrium algorithm 
at each node of the finite element mesh 
OUTPUT: corrected concentration for each ion 
and solid phase profiles 



CHEMICAL DAMAGE f 

Correction of the transport properties 
according to equation (44) 




OUTPUT 

Solid phase profiles 
Ionic concentration profiles 
Adsorbed ion profiles 
Water content profile 
Electrical potential profile 




Figure 2b 
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Figure 4 
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Extended Nemst-Planck model 
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Environment: 
70% Relative Humidity £ f 




Unsaturated soil; 85% relative humidity 
50 mmoVL S04 
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Figure 8 
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Figure 9 
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Electrical coupling and chemical activity 
Electrical coupling without chemical activity 
No electrical coupling nor chemical activity 
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Figure 1 1 
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Figure 12 
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Figure 1 3 



(iCj 50 mmol/L of Na2S04 
in solution 
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Figure 14 



Mesuremencs 
Numerical simulations 
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Figure 15 
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Figure 16 
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Figure 17 
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DETERMINING A FIRST CONCENTRATION 
OF IONS, A WATER CONTENT PROFILE AND 
AN ELECTRICAL POTENTIAL PROFILE USING 
A TRANSPORT ALGORITHM 
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CALCULATING A CORRECTED CONCENTRATION 
OF IONS AND SOLID PHASE PROFILES USING A 
CHEMICAL REACTIONS ALGORITHM 



CALCULATING A CHANGED 
TRANSPORT PROPERTIES PROFILE 
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DETERMINING AN ION CONCENTRATION 
AND SOLID PHASE PROFILES 



FIGURE 18 



(5 



lUUlLLi,..aja-u-jiLuii r ii,. 

MAY, 24. 2001 12:09PM SWABEY OGILVY MTL 514 283 3389 



NO. 0711 r. LULL 



INITIAL DATA 

Material properties: porosity, 

pore solution, 
temperature 
Numerical parameters: number of time steps 

length of time steps, FEM mesh, 
set of tonuosity 
Boudaiy conditions: concentrations 

electrical potential 




TRANSPORT STEP 

Resolution of equations (2), (14) with the finite 
element method v^ith the diffusion coefficients 
corresponding to the tortuosity 
OUTPUT: concentration profile for each ion 
electrical potential profile 



CALCULATE CURRENT 

Calculate the electrical 
current with equation (55) 




OUTPUT 

The numerical current that 
match the measured current 
gives the diffusion 
coefficients for this 
material. 




